Agent-based modeling: a systematic assessment of use cases and requirements for enhancing pharmaceutical research and development productivity.

A crisis continues to brew within the pharmaceutical research and development (R&D) enterprise: productivity continues declining as costs rise, despite ongoing, often dramatic scientific and technical advances. To reverse this trend, we offer various suggestions for both the expansion and broader adoption of modeling and simulation (M&S) methods. We suggest strategies and scenarios intended to enable new M&S use cases that directly engage R&D knowledge generation and build actionable mechanistic insight, thereby opening the door to enhanced productivity. What M&S requirements must be satisfied to access and open the door, and begin reversing the productivity decline? Can current methods and tools fulfill the requirements, or are new methods necessary? We draw on the relevant, recent literature to provide and explore answers. In so doing, we identify essential, key roles for agent-based and other methods. We assemble a list of requirements necessary for M&S to meet the diverse needs distilled from a collection of research, review, and opinion articles. We argue that to realize its full potential, M&S should be actualized within a larger information technology framework--a dynamic knowledge repository--wherein models of various types execute, evolve, and increase in accuracy over time. We offer some details of the issues that must be addressed for such a repository to accrue the capabilities needed to reverse the productivity decline

Propagation of Pericentral Necrosis During Acetaminophen-Induced Liver Injury: Evidence for Early Interhepatocyte Communication and Information Exchange.

Acetaminophen (APAP)-induced liver injury is clinically significant, and APAP overdose in mice often serves as a model for drug-induced liver injury in humans. By specifying that APAP metabolism, reactive metabolite formation, glutathione depletion, and mitigation of mitochondrial damage within individual hepatocytes are functions of intralobular location, an earlier virtual model mechanism provided the first concrete multiattribute explanation for how and why early necrosis occurs close to the central vein (CV). However, two characteristic features could not be simulated consistently: necrosis occurring first adjacent to the CV, and subsequent necrosis occurring primarily adjacent to hepatocytes that have already initiated necrosis. We sought parsimonious model mechanism enhancements that would manage spatiotemporal heterogeneity sufficiently to enable meeting two new target attributes and conducted virtual experiments to explore different ideas for model mechanism improvement at intrahepatocyte and multihepatocyte levels. For the latter, evidence supports intercellular communication via exosomes, gap junctions, and connexin hemichannels playing essential roles in the toxic effects of chemicals, including facilitating or counteracting cell death processes. Logic requiring hepatocytes to obtain current information about whether downstream and lateral neighbors have triggered necrosis enabled virtual hepatocytes to achieve both new target attributes. A virtual hepatocyte that is glutathione-depleted uses that information to determine if it will initiate necrosis. When a less-stressed hepatocyte is flanked by at least two neighbors that have triggered necrosis, it too will initiate necrosis. We hypothesize that the resulting intercellular communication-enabled model mechanism is analogous to the actual explanation for APAP-induced hepatotoxicity at comparable levels of granularity

Dynamics of in silico leukocyte rolling, activation, and adhesion

BACKGROUND: We present a multilevel, agent based, in silico model that represents the dynamics of rolling, activation, and adhesion of individual leukocytes in vitro. Object-oriented software components were designed, verified, plugged together, and then operated in ways that represent the molecular and cellular mechanisms believed responsible for leukocyte rolling and adhesion. The result is an in silico analogue of an experimental in vitro system. The experimentally measured, phenotypic attributes of the analogue were compared and contrasted to those of leukocytes in vitro from three different experimental conditions. RESULTS: The individual in silico dynamics of "rolling" on simulated P-selectin, and separately on simulated VCAM-1, were an acceptable match to individual in vitro distance-time and velocity-time measurements. The analogues are also able to represent the transition from rolling to adhesion on P-selectin and VCAM-1 in the presence of GRO-α chemokine. The individual in silico and in vitro behavioral similarities translated successfully to population level measures. These behavioral similarities were enabled in part by subdividing the functionality of the analogue's surface into 600 independent, "cell"-controlled, equally capable modules of comparable functionality. CONCLUSION: The overlap in phenotypic attributes of our analogue with those of leukocytes in vitro confirm the considerable potential of our model for studying the key events that determine the behavioral outcome of individual leukocytes during rolling, activation, and adhesion. Our results provide an important foundation and framework for future in silico research into plausible causal links between well-documented, subcellular molecular level events and the variety of systemic phenotypic attributes that distinguish normal leukocyte adhesion from abnormal disease-associated adhesion

Characteristics of Oligonucleotide Uptake in Human Keratinocyte Cultures

Oligodeoxyribonucleotides have the potential to interfere selectivity with cellular protein synthesis by sequence-specific hybridization to DNA or RNA molecules. We have, investigated the properties of uptake and intracellular localization of fluorescently labeled oligonucleotides in cultured human keratinocytes using confocal laser scanning microscopy. Unlike many other cell types studied, keratinocytes can internalize oligonucleotides without apparent sequestration in endosomes or cell surface accumulation. Uptake is primarily nuclear and unaltered by sodium azide, monensin, or chloroquin pretreatment. We have verified our results with two different fluorophores, fluorescein and Bodipy, and found similar uptake and distribution patterns in both live and fixed cell populations. Surprisingly, we have found uptake to be heterogeneous within a population, with 15-30% of cells internalizing the oligonucleotides. This percentage is drastically increased to roughly 80% at cell population margins, and after release from M phase arrest. These results on uptake and intracellular localization suggest that keratinocytes may have increased sensitivity as target cells for oligonucleotide based gene regulation strategies

Essential operating principles for tumor spheroid growth

<p>Abstract</p> <p>Background</p> <p>Our objective was to discover in silico axioms that are plausible representations of the operating principles realized during characteristic growth of EMT6/Ro mouse mammary tumor spheroids in culture. To reach that objective we engineered and iteratively falsified an agent-based analogue of EMT6 spheroid growth. EMT6 spheroids display consistent and predictable growth characteristics, implying that individual cell behaviors are tightly controlled and regulated. An approach to understanding how individual cell behaviors contribute to system behaviors is to discover a set of principles that enable abstract agents to exhibit closely analogous behaviors using only information available in an agent's immediate environment. We listed key attributes of EMT6 spheroid growth, which became our behavioral targets. Included were the development of a necrotic core surrounded by quiescent and proliferating cells, and growth data at two distinct levels of nutrient.</p> <p>Results</p> <p>We then created an analogue made up of quasi-autonomous software agents and an abstract environment in which they could operate. The system was designed so that upon execution it could mimic EMT6 cells forming spheroids in culture. Each agent used an identical set of axiomatic operating principles. In sequence, we used the list of targeted attributes to falsify and revise these axioms, until the analogue exhibited behaviors and attributes that were within prespecified ranges of those targeted, thereby achieving a level of validation.</p> <p>Conclusion</p> <p>The finalized analogue required nine axioms. We posit that the validated analogue's operating principles are reasonable representations of those utilized by EMT6/Ro cells during tumor spheroid development.</p

The Effect Of Gender And Knowledge On Students Impressions Of Accountants In The Post-Enron Era

This study examines college students’ impressions of accountants from various sources. It extends previous research on students’ impressions of accountants by analyzing impressions after Enron, gender differences in impressions, and the effect of knowledge on impressions and adding news reports as an impression source. Students rated their impressions of accountants when they selected a major on each of 30 characteristics (such as ethical, versatile, and exciting) from each of four sources, plus a general rating. When qualities were combined into two factors, Professionalism and Personability, accountants were seen as professional, but not particularly personable. In many areas, finance majors had impressions of accountants equal to or higher than those of accounting students, while other majors had significantly lower impressions of accountants. News reports did not significantly influence students’ views of accountants. Knowledge was strongly correlated with impressions. Accounting students had the lowest knowledge of accountants from news reports and movies/TV of any major. Females held more positive impressions of accountants than did males. Implications for various groups, including college and high school accounting education and the accounting profession, are discussed

Computational investigation of epithelial cell dynamic phenotype in vitro

<p>Abstract</p> <p>Background</p> <p>When grown in three-dimensional (3D) cultures, epithelial cells typically form cystic organoids that recapitulate cardinal features of in vivo epithelial structures. Characterizing essential cell actions and their roles, which constitute the system's dynamic phenotype, is critical to gaining deeper insight into the cystogenesis phenomena.</p> <p>Methods</p> <p>Starting with an earlier in silico epithelial analogue (ISEA1) that validated for several Madin-Darby canine kidney (MDCK) epithelial cell culture attributes, we built a revised analogue (ISEA2) to increase overlap between analogue and cell culture traits. Both analogues used agent-based, discrete event methods. A set of axioms determined ISEA behaviors; together, they specified the analogue's operating principles. A new experimentation framework enabled tracking relative axiom use and roles during simulated cystogenesis along with establishment of the consequences of their disruption.</p> <p>Results</p> <p>ISEA2 consistently produced convex cystic structures in a simulated embedded culture. Axiom use measures provided detailed descriptions of the analogue's dynamic phenotype. Dysregulating key cell death and division axioms led to disorganized structures. Adhering to either axiom less than 80% of the time caused ISEA1 to form easily identified morphological changes. ISEA2 was more robust to identical dysregulation. Both dysregulated analogues exhibited characteristics that resembled those associated with an in vitro model of early glandular epithelial cancer.</p> <p>Conclusion</p> <p>We documented the causal chains of events, and their relative roles, responsible for simulated cystogenesis. The results stand as an early hypothesis–a theory–of how individual MDCK cell actions give rise to consistently roundish, cystic organoids.</p

A computational approach to resolve cell level contributions to early glandular epithelial cancer progression

<p>Abstract</p> <p>Background</p> <p>Three-dimensional (3D) embedded cell cultures provide an appropriate physiological environment to reconstruct features of early glandular epithelial cancer. Although these are orders of magnitude simpler than tissues, they too are complex systems that have proven challenging to understand. We used agent-based, discrete event simulation modeling methods to build working hypotheses of mechanisms of epithelial 3D culture phenotype and early cancer progression. Starting with an earlier software analogue, we validated an improved in silico epithelial analogue (ISEA) for cardinal features of a normally developed MDCK cyst. A set of axiomatic operating principles defined simulated cell actions. We explored selective disruption of individual simulated cell actions. New framework features enabled recording detailed measures of ISEA cell activities and morphology.</p> <p>Results</p> <p>Enabled by a small set of cell operating principles, ISEA cells multiplied and self-organized into cyst-like structures that mimicked those of MDCK cells in a 3D embedded cell culture. Selective disruption of "anoikis" or directional cell division caused the ISEA to develop phenotypic features resembling those of in vitro tumor reconstruction models and cancerous tissues in vivo. Disrupting either process, or both, altered cell activity patterns that resulted in morphologically similar outcomes. Increased disruption led to a prolonged presence of intraluminal cells.</p> <p>Conclusions</p> <p>ISEA mechanisms, behaviors, and morphological properties may have biological counterparts. To the extent that in silico-to-in vitro mappings are valid, the results suggest plausible, additional mechanisms of in vitro cancer reconstruction or reversion, and raise potentially significant implications for early cancer diagnosis based on histology. Further ISEA development and use are expected to provide a viable platform to complement in vitro methods for unraveling the mechanistic basis of epithelial morphogenesis and cancer progression.</p

Simulating Properties of In Vitro Epithelial Cell Morphogenesis

How do individual epithelial cells (ECs) organize into multicellular structures? ECs are studied in vitro to help answer that question. Characteristic growth features include stable cyst formation in embedded culture, inverted cyst formation in suspension culture, and lumen formation in overlay culture. Formation of these characteristic structures is believed to be a consequence of an intrinsic program of differentiation and de-differentiation. To help discover how such a program may function, we developed an in silico analogue in which space, events, and time are discretized. Software agents and objects represent cells and components of the environment. “Cells” act independently. The “program” governing their behavior is embedded within each in the form of axioms and an inflexible decisional process. Relationships between the axioms and recognized cell functions are specified. Interactions between “cells” and environment components during simulation give rise to a complex in silico phenotype characterized by context-dependent structures that mimic counterparts observed in four different in vitro culture conditions: a targeted set of in vitro phenotypic attributes was matched by in silico attributes. However, for a particular growth condition, the analogue failed to exhibit behaviors characteristic of functionally polarized ECs. We solved this problem by following an iterative refinement method that improved the first analogue and led to a second: it exhibited characteristic differentiation and growth properties in all simulated growth conditions. It is the first model to simultaneously provide a representation of nonpolarized and structurally polarized cell types, and a mechanism for their interconversion. The second analogue also uses an inflexible axiomatic program. When specific axioms are relaxed, growths strikingly characteristic of cancerous and precancerous lesions are observed. In one case, the simulated cause is aberrant matrix production. Analogue design facilitates gaining deeper insight into such phenomena by making it easy to replace low-resolution components with increasingly detailed and realistic components